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Abstract 

Suppression of T c in bilayer bismuth and thallium cuprates, by substitution 
of Co impurities at Cu sites, are taken for examination. T c suppression data on 
differently doped Bi2212 and T12212 are analysed within the unitary pair-breaking 
formalism due to Abrikosov and Gorkov, by fitting data points to a phenomeno- 
logical relation valid for weak coupling d-wave superconductors. Values of the 
pseudogap magnitude at each doping are thereby estimated within a "fermi-level 
density of states suppression" picture. Pseudogap magnitude from our estimation 
is observed to have a correspondence with a related characteristic temperature T* 
obtained by thermoelectric power measurements. Effects of pseudogap, on the den- 
sity of states, is studied by calculating the susceptibility which shows a broad peak 
at high temperature. This peak feature in susceptibility is indicative of an unusual 
metallic state which could further be explored by systematic other measurements. 
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I. Introduction 

In high-T c layered superconductors (HTLS), studying the suppression of T c by impu- 
rity substitution at the Cu-site of the CuC>2 plane is of considerable importance, because 
it can provide useful information regarding the symmetry of the superconducting order 
parameter, and could also throw light on the recent striking issue of the observation of 
a normal state pseudogap in underdoped materials. In fact, in the recent years, there 
have been growing interest in the studies of suppression of T c in HTLS by intentional 
substitution of both magnetic (Co, Ni) as well as non-magnetic (Zn) impurities at copper 
(Cu) sites. These impurities essentially cause pair breaking |l|, Q leading to the reduc- 
tion of T c while leaving the carrier concentration in the Cu02 plane unaltered [[3], [| . It 
has been shown previously |J that depending on the symmetry of the superconducting 
order parameter, the pair breaking rate or the variation of the experimentally measured 
T c would be different. This is due to the fact that a pairing state of reduced symmetry 
is unstable to impurity scattering to a very high degree, i.e., depression of T c is signifi- 
cantly stronger for a superconductor with d-wave symmetry than for a superconductor 
with pure s-wave symmetry. 

The issue of the symmetry of the superconducting order parameter (OP) in HTLS has 
attracted considerable attention in the recent past and has been discussed widely in the 
high-T c literature, a large fraction of which reported evidence of a d x 2_ y 2 OP symmetry 
[0. 0, and a concensus seems emerging in this direction. But, the striking results relating 
the pseudogap, as reported in NMR || |9], |10|| , optical conductivity [[11], [L2], heat capacity 



[ |l"3|j , transport data |14| and ARPES studies ||15|| , are yet to be understood within an 
unified framework. In other words, there have been no consensus so far regarding the 
origin of the pseudogap and the issue regarding the role played by the pseudogap to affect 
the physical characteristics of HTLS. In this paper, we will primarily be concerned about 
the extraction of the pseudogap and its effects on the in-plane (Cu02 plane) physical 
characteristics. Our basic task is three-fold. First, we analyse the T c -suppression data 
(by impurity substitution) in LM2212 [|[ and T12212 |16[ within the Abrikosov-Gorkov 
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formalism || [T7|, [T8j and investigate the suppression of the in-plane density of states 
(DOS) for variously doped samples. Secondly, we extract the pseudogap magnitude at 



different dopings within a "fermi- level density of states suppression" picture [19|. And 
finally we study the effects of the pseudogap on a measurable physical characteristic of 
HTLS i.e. suceptibility. Note that, T c suppression data, which are obtained by substitu- 
tion of only magnetic impurities (Co) at the Cu-sites, are considered here. An estimation 
of the pseudogap as well as the understanding of its effects on the in-plane physical char- 
acteristics is of importance for the purpose of a consistent theoretical modelling |2(| of 
HTLS as well as for further enlightening of the important issues in the subject. 

The paper is organised as follows. In section-II, we describe the procedural steps in- 
volving the extraction of the pseudogap from the analysis of T c - suppression data. Results 
are presented in section-Ill together with discussions of them. Section-IV comprises of 
a brief summary of the work done and some relevant comments. 

II. Data Analysis: Estimation of the Pseudogap 

The variation of T c as a function of increasing Co impurity concentration has been 
reported elsewhere by B. Bandyopadhyay et al. M. Here, we have fitted the T c - variation 



data to the Abrikosov-Gorkov (AG) equation [17 valid for the weak coupling ci-wave 
superconductors. Within the AG formalism pair breaking rate is inversely proportional 
to the DOS at fermi level denoted by N(Ep). Experimental data are fitted to the 
AG relation using N(Ep) as a fitting parameter and its values are thereby found out 
corresponding to each doping. This DOS at the fermi level N(Ep) shows a progressive 
depletion towards underdoping which could imply the opening up of a gap of growing 
magnitude. Assuming that the depletion in N(Ep) is caused by a pseudogap E g [[15 



we calculate a DOS within the "DOS-suppression picture" using a phenomenological 
form of the quasiparticle energy proposed previously [II], |T^] . Values of the magnitude 



of Eg are then obtained by fitting the calculated DOS at the fermi level (within the 
DOS-suppression picture) to the N(Ep) values found out by data analysis. The E g 



3 



values thus obtained are used to calculate the susceptibility at different doping and 
experimental susceptibility measurement data |HJ are compared. We also do a similar 
analysis, as above, of the T c suppression data on T12212 and report the results here. 

From the analysis of T c reduction data within the AG framework, it is found that 
both Bi2212 and T12212 cuprates bear qualitatively similar characteristic features re- 
garding the pseudogap and other related properties investigated here. It is evident in 
the data analysis that the rate of suppression of T c becomes faster towards underdoping 
as reported earlier by other groups 0, ||]. We find that the magnitude of E g , as obtained 
within the DOS suppression picture, grows as an inverse power of the carrier concentra- 
tion towards underdoped side and is qualitatively similar to a related temperature scale 
T* estimated by thermoelectric power measurements 01]. 

For Bi2212 cuprates the sample composition is given by Bi 2 Sr2Cai_ x Y x (Cui_ y Mj / ) 2 08, 
where M represents cobalt (Co) impurity and y gives the concentration of impurity in per- 
centage. Here x is the amount of yttrium (Y) substituted in place of calcium (Ca) which 
controls the carrier concentration. Details of sample preparation and measurements of 
T c versus y (for Bi2212) are presented earlier 0] and here we consider an analysis of the 
T c suppression data. 

For Tl-cuprate the sample composition is Tl2Ba2Cai_ x Y I .(Cui_ 3/ M J ,) 2 C , 8 where M 
and y again represent Co impurity and its concentration respectively. Values of T c for 



both Bi2212 BJ and T12212 [16] at different impurity levels are estimated by studying 
the temperature variation of the resistivity. It should be mentioned that the carrier 
concentration for the samples used here are changed by substitution of yttrium (Y) at Ca 



sites pi], 021. This procedure of changing carrier concentration does not affect the Cu02 



plane responsible for superconductivity. Whereas a change of the carrier concentration 
by oxygenation procedure might directly affect the Cu02 plane by leaving traces of 
inherent impurities during sample preparation. These impurities are not accounted for 
while studying the T c suppression by impurity substitution and hence, the samples (of 
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varied carrier concentration) prepared by Y-substitution procedure are better suited for 
the purpose. Thus, T c suppression data used here j|, [16] for analysis are expected to be 
more reliable. 



III. Results and Discussions 

In Fig.l we present T c versus y data (symbols) for various values of x (corresponding 
to different doping levels) where Fig. 1(a) includes data for Bi2212 and Fig.l(b) gives that 
for T12212 samples. Clearly, for overdoped situation (small x) the reduction of T c with 
Co-impurity concentration (y) is slower than that for the underdoped situation (large 
x). The T c reduction data are fitted to the AG equation |L7|, |18| given as 



i r 



i) 



.2 2nk B T, 

where T co is the value of T c at y = 0, i.e., without any impurity scattering centre in the 
sample. Here ip[z] represents the digamma function for argument z and T = rii/(jrN(Ep)) 
is the pair breaking scattering rate for unitary scattering. The quantity N(Ep) is the 

o3 

DOS at fermi level (in units of /eV/A ) and n« = ay/abc is the density of impurity 
scatterers per unit volume where a is the number of Cu02 planes per formula unit and 
a, b, c are lattice constants. A justification for using the AG equation for the short 



coherence length cuprate superconductors can be found in references @ and | I3| . In 
Fig. 1(a) and Fig. 1(b) solid and dashed lines are fitted curves using Eq.(l). Matching 
of the fitted lines with the experimental data are quite good, as seen in Fig.l. Notice, 
that the experimental data points (T c versus y) exist only upto a limited value of the 
impurity concentration y = 7 — 8 %. This is because samples of a pure single phase are 
not available f|] beyond this value of impurity concentration, y > 8%. 

In the curve fitting process, N(Ep) is taken as a fitting parameter and its values are 
found out by a least square fitting of T c versus y data. Fitted N(Ep) values corresponding 
to different doping levels (different x) are listed in table-I. It may be noted here that the 
typical r.m.s. deviation in the least square fitting of data is of the order of 1 K which 
reaffirms the visual excellence of data fitting curves presented in Fig.l. 



The ratio |dT c /ck/| is calculated for the fitted curves (T c versus y) of Fig.l. It is seen 
that with the increase of y, the slope |dT c /d?/| increases monotonically for small y and 
then registers a steep rise before falling sharply to zero at some value of y denoted as 
the critical value of the impurity concentration y c . Values of y c are different for different 
doping levels of the samples. For overdoping y c is considerably large and decreases 
progressively but rapidly towards underdoping. Noting the values of y c for each curve, 
one could calculate the corresponding critical values of the pair breaking scattering rate as 
r c = ay c / '(tcJS '(Ep)ab c) which are listed in table-I. The AG relation used here (in Eq.(l)) 
is valid for the weak coupling (i-wave superconductors for which the critical scattering 
rate (at which T c reduces to zero) is evaluated |18| to be r c = 0.88 k B T co . Values of r c , 
calculated by this empirical relationship, are also presented in table-I which are within 
1 % of the values obtained by data analysis. This correspondence might be taken as a 
crosscheck for the data fitting procedure and could also signify the validity of the AG 
relation used for the purpose. 

Next we calculate the DOS using the normal state quasiparticle dispersion given as 
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E* 



el + E g i),f ll \ (2) 



where represents the tight binding band energy for which we consider a realistic 
band structure obtained by a six parameter tight-binding fit [t ^i^2^3^4^5] = [0.131,- 
0.149,0.041,-0.013,-0.014,0.013] eV to the ARPES data on Bi2212 ||. Here E g (k) = 
(Eg/2) |?7k| with E g being the pseudogap magnitude and r] k = cos k x a — cos k y a (a is the 
in-plane lattice constant of a square lattice) is the d x 2_ y 2 symmetry factor associated 



with the pseudogap JT5[. Appearence of E g in the quasiparticle dispersion effectively 
causes suppression of the DOS at the fermi level. In our analysis, we estimate E g at 
a fixed doping level by fitting its value such that the DOS at Fermi energy calculated 
by using relation (2) (within the "DOS-suppression picture") matches N(Ep) found out 
previously from the data analysis (see table-I). A plot of scaled E g (open squares) and T co 
(solid triangles) as a function of carrier concentration (p), for B42212 samples, is given 
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in Fig. 2. Here, E g at different p are scaled to its value at optimal doping p = 0.16, that 
is E g = E g (p) I E g (p = 0.16). T co values obtained in experimental measurements are also 
scaled similarly T co = T co {p) /T co {p = 0.16). Solid line in the main figure, representing 
the locus of T co data, is a guide to the eye and the dashed line is a power law fit 
to the E g data. Notice that the pseudogap magnitude E g shows a sharp rise towards 
underdoping which is seen to vary as a power of inverse doping {1/p). For the case of 
T12212 samples, the data points (not included in the figure) show similar characteristic 
features |24| as in the case of Bi2212. It would be of relevance to point out that in a 



recent experiment, Demsar et al. |25| studied the gap-structure in YBCO single crystals 
employing real-time measurements of the quasiparticle relaxation dynamics. They found 
a T-independent pseudogap to remain dominant for the underdoped samples and even 
persists in the overdoped state with its magnitude being inversely proportional to doping. 

In the inset of Fig.2, we include the variation of T* versus p (solid squares) as esti- 
mated from TEP measurement on Bi2212 where T* is the temperature at which thermo- 
electric power shows a peak and is much higher than the corresponding superconducting 
transition temperature T c . The quantity T* can be thought of as an energy scale related 
to the pseudogap E g , although in some places it has been termed as the temperature 
where the pseudogap opens up J26|. Solid line in the inset is a fit to the data and depicts 
that T* grows towards underdoping as a power of 1/p, although the power is different 
from that involving E g versus p in the main figure. This similarity observed within T* 
and E g emboldens the idea that T* is an energy scale related to E g . It is important 
to note that E g values towards overdoping do not fall to zero, but shows monotonic 
decrease or remain nearly flat implying existence of E g even in the overdoped cuprates. 
This feature is visible in both the main figure as well as in the experimental T* data 
presented in the inset. Similar results of the existence of pseudogap in the overdoped 
region has been reported in a recent experiment on YBCO samples and for the case 
of Bi 2 Sr 2 (Ca,Y)Cu208 as observed in several other experiments |27|, 2l|. This is also 



consistent with earlier experimental observations in Hg-cuprates |29| . 
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With the assumption, that quasiparticles are well defined in the vicinity of the fermi 
surface, susceptibility can be written involving DOS as 

df(EY 



/c 
-I 



N(E) dE (3) 



BE 

where N(E) is the DOS at energy E and f(E) = {er + is the fermi distribution 
function. Susceptibility x becomes temperature dependent through the fermi function 
f(E). Using the quasiparticle dispersion as in Eq.(2), we calculate \ for different carrier 
concentrations (different E g ). A plot of xl V^b as a function of T is given in Fig.3. The 
quantity x/^b nas the dimension of DOS denoted by N(E) in Eq.(3). In each curve 
in the plot, towards low temperature x falls °ff rapidly and tends to zero in the limit 
T —>■ 0. With increasing temperature, x increases slowly, peaks at a temperature (T x ) 
and then decreases slowly towards its high-T saturation. Occurance of the broad peak 
in x & t high T (above the corresponding experimentally measured T co ) is due to the 
presence of E g in the quasiparticle energy and is not observable in case of a normal 
metal. As one goes towards reduced doping level, E g increases, the value of T x gets 
shifted towards higher temperatures and the peak becomes more and more broader. In 



the inset of Fig.3, we include the experimental data [£TJ of x versus T for Bi2212 which 
shows the feature that x has a broad peak above T c and the peak position shifts towards 
higher T with underdoping. This consistency of the earlier experimental data with that 
calculated using the quasiparticle dispersion in Eq. (2) lends some support to the physical 
picture that the existence of E g affects the in-plane (Cu02-plane) charge transport by 
suppressing the quasiparticle DOS at the fermi level. It may be noted here that, within 
a similar DOS suppression picture, NMR data 0, as well as the normal state heat 
capacity and magnetic susceptibility data |JD| on different HTLS materials were analysed 
earlier in terms of a d-wave pseudogap. 

IV. Summary and Comments 

To summarize, we have considered experimental T c -suppression data by impurity 
substitution in Bi2212 and T12212 cuprates, analysed them within Abrikosov-Gorkov 
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formalism of pair breaking by impurity scattering and estimated the pseudogap magnitude 
within the DOS suppression picture. Our analysis shows that the pseudogap energy 
varies as a power of inverse doping and the variation bears close resemblance to that 
of T* found by thermoelectric power measurement on Bi2212. Presence of E g affects 
the susceptibility x by inflicting the occurance of a broad peak above the corresponding 
experimental T co indicating an unusual metallic state. A systematic study of x an d other 
properties connected to the pseudogap, for differently doped samples of various layered 
cuprates could be useful to further explore the relevance of the DOS suppression (by 
pseudogap) viewpoint in layered cuprates. 
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Table-I 



Material x p N(Ep) T c values (in meV) 

o3 

(/eV/A ) From curve fitting From the relation 

r c =0.88 k B T co 



Bi2212 


0.0 


0.225 


0.0805 


5.90 


5.86 




0.2 


0.160 


0.0435 


6.77 


6.73 




0.3 


0.135 


0.0250 


7.09 


6.98 




0.4 


0.110 


0.0191 


6.09 


6.03 


T12212 


0.2 


0.102 


0.0547 


5.94 


5.89 




0.25 


0.085 


0.0466 


4.19 


4.15 
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Figure captions: 



Fig.l. Superconducting transition temperature T c is plotted as a function of substituted 
Co impurity concentration for different values of yttrium content x as shown in the 
figure. Figure (a) corresponds to Bi2212 and (b) corresponds to T12212 samples. 
In both the panels, the symbols represent experimentally measured T c values and 
lines (solid and dashes) represent fitted curves using the AG equation (1) as written 
in the text. 

Fig.2. Phase diagram of bilayer bismuth cuprate. Solid triangles are scaled T co values 
corresponding to different carrier concentration (p) as measured experimentally. 
Solid line is a parabolic fit to these experimental data and serves as a guide to the 
eye. Open squares are scaled values of the pseudogap magnitude (E s ) as estimated 
within the DOS suppression picture. Dashed line is a power law fit to E g data. 
[Inset: Solid squares represent values of T*, as estimated by thermoelectric power 
measurement, as a function of p. Solid line is a power law fit to the data.] 

Fig.3. Susceptibility for Bi-cuprate, as calculated using Eq.(3) of the text, is plotted as 
a function of temperature for three different values of x mentioned in the figure. 
The temperature Tp for each curve, at which susceptibility has a broad peak, is 
marked by an arrow and the corresponding value is written. [Inset: Susceptibility 
data for Bi2212 as a function of temperature corresponding to three different values 
of yttrium (Y) concentration as noted in the figure. These data are taken from 
Ref.jnj.] 
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